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Abstract. We have assessed the influence of the stellar iron content on the Cepheid Period-Luminosity (PL) relation by relating 
the V band residuals from the Freedman et al <2001> PL relation to [Fe/H] for 37 Galactic and Magellanic Clouds Cepheids. 
The iron abundances were measured from FEROS and UVES high-resolution and high-signal to noise optical spectra. Our 
data indicate that the stars become fainter as metallicity increases, until a plateau or turnover point is reached at about solar 
metallicity. Our data are incompatible with both no dependence of the PL relation on iron abundance, and with the linearly 
decreasing behavior often found in the literature (e.g. Kennicutt et al ll998l Sakai et al l2004ii . On the other hand, non-linear 
theoretical models of Fiorentino et al <2002t provide a fairly good description of the data. 
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1. Introduction 

. Ever since the work of Edwin Hubble, the Cepheid Period- 

' Luminosity (PL) relation is a fundamental tool in determining 
Galactic and extragalactic distances. In spite of its paramount 
importance, to this day we still lack firm theoretical and em- 

i pirical assessment on whether or not chemical composition has 

' any significant influence on it. 

Theoretical pulsational models by different groups lead to 
markedly different results. On the one side computations based 
on linear models (e.g. Chiosi et al l 19921 Sandage et al 1999 
Baraffe & Alibert .2001 J suggest a mild dependence of the 
PL relation on chemical composition. The predicted change at 
log(P) = 1 is less than 0. 1 mag at all wavelengths between the 
V and K bands for a change in metallicity from Z - 0.004, 
representative of the Small Magellanic Cloud (SMC), to 0.02, 
typical of Galactic Cepheids at the Solar circle. This result is 
challenged by the outcome of the non-linear convective mod- 
els (e.g. Bono et aI II999l Caputo et aI I2000> . which find that 
both the slope and the zeropoint of the Period-Luminosity re- 
lation depend significantly on the adopted chemical composi- 
tion. Again for log(P) = 1 and the same variation in metallicity 
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as above, they predict a change as large as 0.4 magnitudes in 
V, 0.3 magnitudes in / and 0.2 magnitudes in K. Moreover, 
the change is such that metal-rich Cepheids are fainter than 
metal-poor ones, again at variance with the results of Baraffe & 
Alibert <2001l l. Recent calculations by Fiorentino et al (20(7^. 
also based on non-linear models, indicate that the PL relation 
also depends on the helium abundance. 

Observationally, indirect measurements in external galax- 
ies from secondary abundance indicators tend to find that 
metal-rich Cepheids are brighter than metal-poor ones, al- 
beit with disappointingly large range of quoted values any- 
where between (Udalski et al Ciardullo et al EUU^ 
and -0.9 mag/dex (Gould lT994l l. For lack of better evidence, 
a marginally significant dependence of -0.2 + 0.2 mag/dex 
was adopted in the final paper of the hst Key Project on Hq 
(Freedman et al '2001'). This was largely based on the results 
of Kennicutt et al (1998 1 from oxygen abundances of Hii re- 
gions in two Cepheid fields in MlOl (-0.24 + 0.16 mag/dex) 
and in 10 Cepheid galaxies with Tip of the Red Giant Branch 
distances (-0.12 + 0.08 mag/dex). This latter method was re- 
cently applied by Sakai et al (2004i to a larger sample of 17 
galaxies yielding -0.24 + 0.05 mag/dex. Storm et al (20041 
assumed a metallicity difference between the Milky Way and 
the Small Magellanic Cloud of 0.7 dex to derive a slope of 
-0.21 + 0.19 mag/dex from 5 SMC and 34 Galactic Cepheids 
for which they measured the distances with a Baade-Wesselink 
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technique. Recently Groenewegen et al ( I2004I I have assembled 
from the literature and homogenized a large sample of Cepheid 
stars with both distance and metallicity determinations to de- 
rive a metallicity effect of -0.27 + 0.08 mag/dex in the zero 
point in VIWK bands. 

Compared to these previous studies, the novelty of our ap- 
proach consists in measuring directly the chemical composi- 
tion of Cepheid stars with know a distance, without relying on 
proxies such as oxygen nebular abundances derived from spec- 
tra of H II regions at the same galactocentric distance of the 
Cepheid field (e.g. Kennicutt et al 1998 1, or secondary distance 
indicators like the Tip of the Red Giant Branch (e.g. Sakai et al 
120041 . 

The paper is organized as follows. The observations are 
presented in section |2 together with a brief description of the 
[Fe/H] measurements. In section |3] the dependence of the PL 
relation on [Fe/H] is derived. Finally, section |4] contains the 
discussion and the conclusions. 

2. Observations and data reduction 

In order to achieve an extended coverage in metallicity, we 
have observed Cepheid stars in three environments known to 
have significantly different mean chemical compositions: the 
Solar Neighborhood, the Large Magellanic Cloud (LMC) and 
the Small Magellanic Cloud. 

The spectra of the 13 Galactic stars were obtained with the 
FEROS instrument (Pritchai'd 2004') at the ESQ 1.5m tele- 
scope on Cerro La Silla. The spectral resolution is 48,000 
and the signal-to-noise ratio is about 70 to 150, depending on 
the brightness of the target and on the spectral range. The 12 
Cepheids in the LMC and the 12 Cepheids in the SMC were 
observed with the UVES spectrograph (Kaufer et al l2004n at 
the VLT-Kueyen telescope on Cerro Paranal. For them the res- 
olution is about 40,000 and the signal-to-noise ratio about 50 
to 70. The 2-D raw spectra were run through the respective in- 
strument pipelines, yielding to 1-D extracted, wavelength cal- 
ibrated and rectified spectra. The normalization of the contin- 
uum was refined with the IRAF task continuum. The 1-D spec- 
tra were corrected for heliocentric velocity using the rvcorr and 
dopcor IRAF task. This latter task was also used to apply the 
radial velocity correction, which was derived from 20 Fel, Fell 
e Mgl lines. 

Selected characteristics of the programme stars are listed in 
Table 1. 

2.1. Metallicity determination 

The complete description of the method we have used to derive 
the metallicity for our programme stars will be presented in a 
forthcoming paper (Mottini et al, in preparation). Here we just 
summarize the salient points. 

The iron content of the Cepheids was derived from the 
equivalent widths (EWs) of 150-200 Fe I and 10-15 Fe II 

' see also 
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unblended lines. The EWs themselves were measured semi- 
interactively with a software developed by one of us (PF, fit- 
line). For about 15% of the lines it was necessary to use 
the splot task in IRAF, instead, because the Gaussian shape 
adopted by fitline could not satisfactorily reproduce the ob- 
served profile (e.g. very broad or asymmetric lines). The equiv- 
alent widths used to determine the iron content range approxi- 
mately from 5 to 150 mA, well sampling the linear part of the 
curve of growth. Full details on the selection of the iron lines 
and their physical properties (oscillator strengths, etc.) will be 
given in Mottini et al (in preparation). 

Following Kovtyukh & Gorlova (2000 1, we have deter- 
mined the stellar effective temperature from 32 line-depths ra- 
tios. Gravity and microturbulent velocity were constrained by 
imposing the ionization balance and by minimizing the slope 
of log(e(Fe)) vs EW, respectively. The LTE stellar model at- 
mospheres by Kurucz (1993) and the WIDTH9 code (Kurucz 
IT993> were used throughout the analysis. The detailed descrip- 
tion of this procedure is beyond the scope of this Letter and 
will also be presented in Mottini et al (in preparation). 

The mean value [Fe/H] is about solar for our Galactic sam- 
ple, 0.4 for the LMC sample and ~ -0.7 for the SMC one, 

with an rms of about 0.15 dex. 

2.2. Periods, photometry and distances 

We have adopted the periods, distance moduli and V-band pho- 
tometry of the Galactic Cepheids as listed in Table 3 of Storm 
et al (2004 l. Two of our programme stars, ^ Gem and y6 Dor, are 
not included in that list and for them we have used the values 
from Groenewegen et al (2004 Table 3). 

The periods and V-band photometry for the Magellanic 
Cloud Cepheids were taken from Laney & Stobie (1994i. The 
distance modulus of the barycenter of the LMC is assumed to 
be 18.50, for consistency with the PL relation of Freedman 
et al (2001 see below). The SMC is considered 0.44 magni- 
tudes more distant (e.g. Cioni et al 2000 1. Depth and projec- 
tion effects in the Magellanic Clouds were corrected for us- 
ing the position angle and inclination of each galaxy as deter- 
mined by van der Marel & Cioni (.2001. LMC) and Caldwell & 
Laney ( fT99T] SMC). 

The periods of our programme stars range between about 
5 and 65 days, thus populating the linear part of the PL re- 
lation, the one useful for distance determinations (e.g. Bono 
et al '1999V All of stars are bona fide fundamental mode pul- 
sators. In particular, the two stars with periods shorter than 8 
days (V Cen and HV 6093) follow the fundamental mode PL 
relation of Freedman et al (2001), with deviations from it which 
are consistent with its intrinsic scatter (0.2 and less than -0.05 
magnitudes, respectively). 

3. The PL relation vs [Fe/H] 

In Figure we plot S(Mv), the V-band residuals of our pro- 
gramme stars from the standard PL relation of Freedman et 
al (I2001> . as a function of the iron abundance we have derived 
from the FEROS and UVES spectra (a positive S{Mv) means 
fainter than the mean relation). The Freedman et al (.200 Ij PL 
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relation was derived for the LMC as a whole ([Fe/H] =^ -0.4) 
and 6{Mv) is the correction to be applied to it as a function of 
metallicity. 

The data in Figure ^ are binned in metallicity to reflect the 
typical uncertainty on our determination of [Fe/H], marked by 
the horizontal errorbars. The median value of 6{Mv) in each 
metallicity bin is plotted as filled squares, with the vertical er- 
rorbars representing its associated errors. The rms on 5{Mv) in 
each bin is of the order of 0.3 mag, corresponding to the intrin- 
sic width of the instability strip. 

As it can be seen in Figure ^ our data indicate that 
6(Mv) increases with [Fe/H] up to about solar metallicity, 
i.e. Cepheids become fainter as metallicity increases, where it 
shows a turnover or a flattening. Regrettably, the rather large 
errorbars on 6(Mv) do not allow us to distinguish between 
these two possibilities. On the other hand, the data are only 
marginally consistent with 6(Mv) still rising at metallicities 
higher than solar: a linear extrapolation of the trend defined by 
the three points at lower metallicity would be 1 .6cr away from 
the measured value at [Fe/H] = 0.15. 

In order to enhance the statistics in the highest metallic- 
ity bin, we have included in our sample 3 stars, SZ Aql, WZ 
Sgr and KQ See, with published metallicities (Andrievsky et 
al I2002al l2002bl [Fe/H] ^ 0.15) and distances, which we 
have derived by combining the radii measured by Laney & 
Stobie ([T995:( with K-hand photometry (Laney & Stobie 1994) 
and the surface brightness calibration of Groenewegen (2004i. 
The result is shown as an open square in Figure^ With this ad- 
dition, the fact that 6(Mv) keeps rising up to super-solar metal- 
licity can be excluded at the 3cr level. 

3.1. Comparison with previous results 

Let us now use a technique to compare our data to selected 
previous results: 

- null hypothesis, i.e. no dependence of SiMy) on [Fe/H]: 
the test returns a value of 13.5. For 3 degrees of free- 
dom, this means that the hypothesis can be discarded at the 
99.6% level. The data can be reconciled with the null hy- 
pothesis if a distance modulus of 18.2-18.3 is used for the 
LMC (x^ ^3.1 for 3 degrees of freedom). Such a short 
distance, however, seems nowadays to be rather implausi- 
ble (e.g., RomanieUo et al EUUni Walker OUU^ Alcock et al 
L2004' Borissova et al"200Tl 

- monotonically decreasing 6iMv) (e.g. Kennicutt et al 
[19981 Sa kai et al l2004l Storm et al l2U04l Groenewegen et 
al "2004*1: notably, this sort of behavior was used in the final 
paper of the hst Key Project on Ho (Freedman et al 200Ll. 
As an illustration, we report in Figure^the classical results 
of Kennicutt et al (WW 6(Mv) = (-0.24 + 0.16) x [Fe/H], 
open circles and solid line). The;^'^ test on it returns a value 
of 28, meaning that it is incompatible with our data with a 
confidence higher than 99.95%. 

- theoretical models by Fiorentino et al 12002> : fliese 
are the outcome of non-linear pulsation computations and 
predict a non-monotonic behavior of 6{Mv) with [Fe/H]. 
Following the prescriptions in that paper, we have com- 



puted the value of 6(Mv) in each metallicity bin by com- 
bining the theoretical Period-Luminosity-Color (PLC) re- 
lations as a function of metal content (see their Table 4) 
with the observed periods and V - K colors of our stars. 
The use of the PLC relations supplies individual absolute 
magnitudes, thus avoiding deceptive uncertainties due to 
the intrinsic width of the instability strip. In the left panel 
of Figure [2 we report as a dashed line the results for the 
models with a helium-to-iron enrichment of AY/AZ = 2.5. 
The x^ value for this model is 4.2 which, for 3 degrees of 
freedom, indicates a fairly good agreement with the data. 
The other two models presented in Fiorentino et al ( 2002>. 
AY/AZ = 3 and 3.5, also result in similarly acceptable val- 
ues ofx^- 

4. Discussion and conclusions 

We have assessed the influence of the stellar iron content on 
the Cepheid Period-Luminosity (PL) relation by relating the V 
band residuals from the Freedman et al (1200 U PL relation to 
[Fe/H] for 37 Galactic and Magellanic Clouds Cepheids. The 
iron abundance was measured from FEROS and UVES high- 
resolution and high-signal to noise optical spectra. The novelty 
of our approach consists in measuring directly the chemical 
composition of Cepheid stars with a known distance, without 
relying on proxies such as oxygen nebular abundances derived 
from spectra of H II regions at the same galactocentric distance 
of the Cepheid field (e.g. Kennicutt et al 1998 ) , or secondary 
distance indicators Uke the Tip of the Red Giant Branch (e.g. 
Sakai et al'2004). Here we have concentrated on the influence 
of the iron content, but the very same analysis can be applied 
to any other chemical species (Mottini et al, in preparation). 

Our results are summarized in Figure ^ (filled squares), 
together with the empirical results of Kennicutt et al C1998.) 
in two Cepheid fields in MlOl (open circles and solid line) 
and the theoretical predictions by Fiorentino et al (2002) from 
non-linear pulsational models (dashed line). The data indicate 
that 6(Mv), the correction to a metal-independent PL rela- 
tion, increases as the iron content increases, i.e. the stars be- 
come fainter as the metallicity increases, until a flattening or 
a turnover is reached at about solar metallicity. The possibil- 
ity that the increasing trend continues at higher metallicities, 
while it cannot be ruled out completely by the current data, is 
disfavored at the 1 .6cr level. 

Applying ax^ technique, the null hypothesis, i.e. no depen- 
dence of the PL relation on the iron content, can be excluded at 
the 99.6% level. Also, empirical monotonically decreasing lin- 
ear relations (e.g. Kennicutt et ari998 Sakai et al'2004' Storm 
et al .2004. Groenewegen et al 2 004) are incompatible with 
our data with a confidence level higher than 99.95%. A better 
agreement is found with the theoretical models of Fiorentino et 
al (2002 ) , which do predict a non-monotonic behavior with a 
turnover at about solar metallicity. The x^ value is 4.2 (3 de- 
grees of freedom) for their model with AY/AZ - 2.5, plotted as 
a dashed line in Figure ^ A similarly good agreement is found 
for their models with AY/AZ = 3 and 3.5. 

It is apparent from an inspection of Figure [Qthat, because 
of the limited number of stars, the errorbars on the data are 
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[Fe/H] [Fe/H] 

Fig. 1. The V-band residuals compared to the Freedman et al f2001^ PL relation are plotted against the iron content measured from 
FEROS and UVES spectra. Left panel: results for the individual stars (grey symbols: Galaxy as circles, LMC as star symbols, 
SMC as squares) and median value in each metallicity bin (filled black squares) with its associated errorbar. Right panel: Filled 
squares represent the median value in each metallicity bin, with is associated errorbar The open square includes the three stars 
with metallicity from Andrievsky et al J2002al i2002b.) and radius from Laney & Stobie ( 1 995 1. The metallicity dependence as 
inferred by Kennicutt et al il998 1 from two Cepheid fields in MlOl (open circles) is shown as a full line. The dashed line shows 
the theoretical predictions by Fiorentino et al (2002) for a helium-to-metal enrichment of AF/AZ = 2.5 and the observed period 
and V - K distributions of the programme stars. 



quite large. However, the error on 5{Mv) scales as the square 
root of the number of stars in each bin. It is, then, just a matter 
of gathering a larger sample of stars and apply the method we 
have outlined here to characterize in detail the dependence of 
the Cepheid PL relation on the stellar chemical composition. 
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Table 1. Selected properties of the programme Cepheid 
stars. Iron abundances were measured from our FEROS and 
UVES spectra (see section 12. H . Periods, distance moduli 
and V-band photometry of the Galactic Cepheids are taken 
from Storm et al J2004> . except f Gem and p Dor from 
Groenewegen et al ( 2004 1. The periods and V-band photometry 
for the Magellanic Cloud Cepheids were taken from Laney & 
Stobie (0994 1. The distance modulus of the barycenter of the 
LMC is assumed to be 18.50 (see text). The SMC is considered 
0.44 magnitudes more distant (e.g. Cioni et al l2000> . Depth and 
projection effects in the Magellanic Clouds were corrected for 
using the position angle and inclination of each galaxy as de- 
termined by van der Marel & Cioni (.2001. LMC) and Caldwell 
& Laney ( fT9^ SMC). 



Name 


log(P) 


Mv 


FFe/Hl 




Galaxy 




V Cen 


0.740 


-3.295 


+0 04 


KN Pen 


1.532 


-6.328 


+0.17 


V VV \^C11 


1.177 


-4.037 


n CO 

u.uz. 


S Nnr 


0.989 


-4.101 


4-0 09 


T Mon 


1.432 


-5.372 




UNor 


1.102 


-4.415 


+0.07 


AQ Pup 


1.479 


-5.513 


-0.07 


VZ Pup 


1.365 


-5.009 


-0.17 


RS Pup 


1.617 


-6.015 


+0.09 


RZ Vel 


1.310 


-5.042 


-0.19 


1 Car 


1.551 


-5.821 


00 


hpta l~)nr 


0.993 


-3.920 


-0.14 




1.007 


-3.897 


-0.19 




LMC 






n V i 1/ i J 


1.382 


-5.037 


U.J / 


n V iuz.j 


1.425 


-4.994 


94. 


n V iz,H-jz. 


0.941 


-3.899 


f) 9Q 


uv 1 7700 


0.911 


-3.627 


n 9Q 


HV 2260 


1.112 


-4.067 


-0.34 


HV 2294 


1.563 


-6.050 


-0.46 


HV 2352 


1.134 


-4.656 


-0.43 


HV 2369 


1.684 


-6.215 


-0.64 


HV 997 


1.119 


-4.308 


-0.14 


HV 6093 


0.680 


-3.338 


-0.49 


HV 2580 


1.228 


-4.833 


-0.18 


HV 2733 


0.941 


-4.164 


-0.26 




SMC 






HV 11211 


1.330 


-5.296 


-0.67 


HV 1365 


1.094 


-4.147 


-0.67 


HV 1954 


1.223 


-5.325 


-0.67 


HV 2064 


1.527 


-5.440 


-0.55 


HV 2195 


1.621 


-5.874 


-0.61 


HV 2209 


1.355 


-5.519 


-0.56 


HV 817 


1.277 


-5.348 


-0.74 


HV 823 


1.504 


-5.336 


-0.65 


HV 824 


1.818 


-6.669 


-0.71 


HV 837 


1.631 


-5.842 


-0.69 


HV 847 


1.433 


-5.282 


-0.67 


HV 865 


1.523 


-6.015 


-0.79 



